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Abstract. Substructure in galaxy clusters can be quan- 
tified with the robust A statistics (Dressier and Shectman 
1988) which uses velocity kinematics and sky projected 
positions. We test its sensitivity using dissipationless nu- 
merical simulations of cluster formation. As in recent ob- 
servations, about 30% of the simulated clusters show sub- 
structure, but the exact percentage depends on the cho- 
sen limit for defining substructure, and a better discrimi- 
nator is the distribution function of the A statistics. The 
Dressler-Shectman statistics correlate well with other sub- 
cluster indicators, but with large scatter due to its sensi- 
tivity to small infalling groups and projection effects. 
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1. Introduction 

During the last two decades, substructure was detected 
in a significant fraction of groups and clusters of galaxies, 
both as double or secondary maxima in the sky-projected 



galaxy distribution of clusters (cf. Baier 1979 and Geller 



& Beers 1982), and as deviations of spherical symmetry in 



X-ray contour maps (Forman et al. 1981). However, quan- 



tifying substructure in clusters is a non-trivial problem. 
Analyzing substructure indicators with numerical simula- 
tions, West et al. (1988) concluded that many subclumps 
may be pure chance projections, and that the remain- 
ing abundance of substructure can barely discriminate be- 
tween a wide range of cosmological scenarios. Projection 
effects are less pronounced in X-rays, and recent X-ray 
data have provided more reliable information on the abun- 
dance of substructure (cf. Mohr et al. 1993| and Buote & 
Xu 1997). 

In the optical spectral range, one gets better results by 
including velocity information for a large sample of clus- 
ter galaxies. Substructure was identified using subgroup 
velocity statistics in clusters ( Dressier fc Shectman 1988| ), 
finding velocity offsets of cD galaxies with respect to 



the major clust er (Bird 1994), employ ing the hierarchical 
tree algorithm (Serna fc Gerbal 1996 ) and wavelet analy- 
sis ( |Girardi et al. 1997 ). Recently, Dressier & Shectman's 
statistics were used by Zabludoff & Mulchaey (1998) for 
6 poor groups, and by Solanes et al. (1998) for 67 rich 
clusters from the ENACS survey. These different analyses 
agreed in (30 — 40)% of clusters which showed statistically 
significant substructure. But it appears that this amount 
depends on the analysis method ( Pinkncy et al. 1996 ) and 
on the imposed reliability criterion. 

The existence of substructure in clusters suggests that 
the clusters are young objects since loosely bound sub- 
groups can survive only a few cluster crossing times, i.e., 
shorter than the Hubble time. In the framework of the 
hierarchical structure formation scenario, clusters grow 
mainly by mergers of smaller objects and by accretion; 
substructure point to recent major mergers. Starting from 
these ideas and the early termination of growth of density 
perturbations in a low density universe, the abundance 
of substructure was suggested as an effective measure of 
the mean matter density in the universe (Richstone, Loeb 
& Turner 1992, Bartelmann , Ehlers & Schneider 1993, 
Kauffman & White 1993, and Lacey & Cole 1994). The 
interpretation of subclusters as recent mergers might ex- 
plain also some of the properties of cD cluster galaxies, 
such as their orientation with respect to the environment 



(West 1994) and the peculiar velocity distribution in rich 



clusters (Mcrritt 1985) 



The velocity kinematics as a signature of substructure 
traces the galaxy distribution, therefore we employed dis- 
sipationless simulations in a large cosmological environ- 
ment. The cosmological models take COBE-normalized 
perturbation spectra in four cosmological scenarios. Ear- 
lier theoretical studies concentrated on the analysis of 
the density contrast of clusters and the density profile 
flJiug et al. 1994 [Thomas et al. 1997|). There, relatively 
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small differences in the cluster properties are found when 
comparing cluster profiles at the same overdensity. This 
concerns mainly the properties of relaxed clusters. Here 
we study in particular clusters which represent deviations 
from an equilibrium. For quantifying substructure in X- 
ray profiles of clusters, hydrodynamic simulations of clus- 
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ter formation have to be employed, and early attempts at 



this have already yielded first promising results (Cen 1997 
and Valdarnini, Ghizzardi & Bonometto 1999). 

The outline of the paper is as follows. First we ana- 
lyze the method for identification of substructure in galaxy 
clusters selected from numerical simulations. Then we ap- 
ply the algorithm to different cosmological models. In Sec- 
tion 4 we compare the substructure measure with other 
methods and with observational data. We conclude with 
a discussion of our results. 

2. Method 

The Dressier & Shectman (1998) statistics evaluate the 
velocity kinematics of galaxy groups identified in sky pro- 
jected clusters. To be specific, one takes a number of neigh- 
bours N nn from each galaxy in the projection, determines 
the mean velocity vi OC ai and velocity dispersion a Vi i oca \ of 
the subsample, and compares this with the mean velocity 
v and velocity dispersion o~ v of the whole group, 



N, 



$i = [("local - v) 2 + (o-„,local _ °~v) 2 ] 



(1) 



A measure of the amount of dumpiness in the cluster 
is the sum of the individual positive Si over all cluster 
galaxies N, 



N 



(2) 



which is called the delta-deviation. The A-deviation is 
large for groups with kinematically distinct subgroups. 
Fig. |l| illustrates the statistics for a simulated typical dark 
matter cluster selected with a linking length of 0.2 times 
the mean interparticle spacing. Around each point a circle 
is plotted with radius proportional to exp(dj) . The cluster 
has a pronounced centre indicated by the central parti- 
cles in the figure decorated with the small circles. To the 
right and above from the centre there are some subclumps 
which do not represent especially important subgroups 
with decoupled kinematics, but most probably they are 
small satellites slowly falling onto the group centre. To 
the left, there are some particles with large Si that con- 
tribute mostly to the cumulative A-deviation. The cluster 
in Fig. |l| represents an example for marginally reliable sub- 
structure, we base our later statistical analysis on clusters 
with more pronounced substructure. An obvious advan- 
tage of the method is that no a priori selections or as- 
sumptions about the positions of subclumps have to be 
imposed. 

In order to test and calibrate the statistics, we have 
compared five clusters of masses of about lO^/i^M© 
(1000 particles), 5 ■ 10 14 /i _1 Mq (500 particles), and 2 • 
10 14 /i _1 M Q (200 particles) in a standard CDM simula- 
tion. For each mass we study clusters both with and with- 
out obvious substructure as decided by eye from the 3- 
dimcnsional distribution. 
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Fig. 1. Bubble plot for a particle group containing sub- 
structure. 

The A-deviation statistics depend on the number of 
cluster galaxies with measured redshift N that can be in- 
cluded in the analysis, and on the number of neighbours 
in tested subgroups N nn . If one assumes a random distri- 
bution of velocities one expects values A oc N (Dressier 
& Shectman 1 988), and for galaxy clusters with substruc- 
ture, A > N (Pinkney et al. 1996). Therefore we normal- 
ize A with the total number of galaxies investigated, and 
in Fig. U we plotted A/N against the number of neigh- 
bours N nn in subgroups. The selected groups with sub- 
structures (solid lines) show increasing values if the clus- 
ters have a large mass. Only the poor clusters lead to a 
poor discrimination of substructure. The curves suggest 
the use of a rather high value for the number of neigh- 
bours, while Pinkney et al. (1996) propose N nn = N 1 ^ 2 . 
The maximum for low mass objects (more than 200 parti- 
cles) lies at about 25 neighbours which we therefore took 
as a reasonable value for the number of nearest neigh- 
bours in the following analysis. We took a fixed number 
in order to be independent of the sampling rate of veloci- 
ties. A similar plot of A/N over the number of objects in 
clusters randomly selected demonstrates a very stable dis- 
crimination of clusters with substructure from those that 
are in equilibrium, and A/N remains almost constant for 
values N > 100. Ideally this would be the number of re- 
quired redshifts to be measured for analyzing clusters for 
substructure. 



3. Application to Cosmological Simulations 

We apply the substructure statistics to a set of cos- 
mological models that are currently under investigation. 
Three of the models are COBE normalized according to 
the prescription of Bunn & White (1997) where we as- 
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Fig. 2. Variation of A/N with the number of neighbours 
AT nn . The thickness of lines is proportional to the mass. 

Table 1. Cosmological scenarios and percentage of sub- 
structure. The box size L is given in /i _1 Mpc, and the 
particle mass m p in units of 10 fa M©. 





«o 


h 


o"8 


L 


m p 


P(A/N > 1.4) 


SCDM1 


1.0 


0.5 


1.18 


200 


11 


33% 


SCDM2 


1.0 


0.5 


0.53 


200 


11 


43% 


ACDM 


0.3 


0.7 


1.00 


280 


9 


27% 


OCDM 


0.5 


0.7 


0.96 


280 


15 


30% 



sume pure adiabatic perturbations and a baryon con- 
tent of tt b h 2 = 1.3 • 10~ 3 . The standard CDM model 
(SCDM1) is taken as a reference model despite its en- 
hanced power at cluster scales. In particular it produces 
too high a cluster abundance. As more realistic alterna- 
tives we take an OCDM model with Q a = 0.5 and a ACDM 
model with f2 = 0.3 and a cosmological constant to pro- 
vide spatial flatness (i.e. FIao — 0.7). Both models are 
promising since they lead nearly to the observed cluster 
abundance (Eke et al. 1996), and they represent models 
which reproduce the observed superclustering of galax- 
ies ( Doroshkevich et al. 1999). In addition to these three 
models we have used the output of the SCDM1 model 
at a redshift z = 1.0 to ensure the correct cluster abun- 
dance (SCDM2). The model parameters are summarized 
m Table § In particular it is shown that the simulations 
are running in boxes that are identical in physical size 
L, thereby making small differences in the mass resolu- 
tion TOp despite the different density parameters. All mod- 
els were run using the adaptive P 3 M code of Couchman 
(1991) with 128 3 particles. An analysis of the virialization 
of clusters in these simulations can be found in a previous 
paper (Knebe & Miiller 1999). 
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Fig. 3. Cumulative probability distribution of A for dif- 
ferent cosmological models. 



As explained in the introduction, we expect differences 
in the abundance of substructure for different cosmologies, 
and therefore in the probability distribution of the delta- 
deviation. In a low-density universe (open or flat), struc- 
ture formation ceases at earlier times compared to the 
SCDM models. This means that clusters in low-density 
universes should show less substructure since they formed 
earlier and therefore had more time to virialize. We al- 
ways obtained strongly varying delta-deviations, there- 
fore, we only analyze its abundance distribution. In Fig. |^ 
we show the cumulative probability distributions of the 
delta-deviations A/N of friends-of- friends clusters selected 
with dimensionless linking length (in terms of the mean in- 
terparticle separation) 11 = 0.2, 0.17, and 0.16 for SCDM, 
OCDM, and ACDM, respectively. We show these distri- 
butions for a fixed number density of simulated clusters 
n c i = 10- 5 fa 3 Mpc~ 3 . 

The distributions show differences which underline the 
cited expectation. In particular, both SCDM simulations 
have higher probabilities of substructure than the low den- 
sity models. Especially the SCDM2 model shows a large 
number of clusters with A/N > 1.4. On the other hand, 
the low density models ACDM and OCDM have a simi- 
lar distribution of the delta-deviation, with ACDM lying 
below the OCDM model in agreement with the expecta- 
tion due to the lower density parameter. To quantify the 
results, the probability of substructure is given in Table [l] 
as the percentage of clusters with A/N > 1.4. Wc find 
that about 30% of the investigated objects show a sig- 
nificant level of substructure. This coincides with the ob- 
served probability of substructure in recent studies using 
the delta-deviation (Dressier & Shectman 1988, Zablud- 
off , 



Mulchaey 1998, Solanes, Salvado-Sol'e & Honz'alez- 
Casado 1999). To check the significance, a Kolmogorov- 
Smirnov test was employed that is summarized in Ta- 
ble 0. The first column gives the maximum distance of 
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Table 2. Kolmogorov-Smirnov (KS) test for the cumu- 
lative distributions. First column for each redshift gives 
maximum distance D, and second column the significance 
level. 



Simulations 


difference 


reliability 


SCDM1 vs. SCDM2 


0.14 


89 % 


SCDM1 vs. ACDM 


0.24 


96 % 


SCDM1 vs. OCDM 


0.18 


95 % 


SCDM2 vs. ACDM 


0.27 


99 % 


SCDM2 vs. OCDM 


0.17 


97 % 


ACDM vs. OCDM 


0.11 


48 % 



0.5 1.0 1.5 2.0 2.5 

A / M 



the compared distributions, whereas the second column is 
the probability (in percent) that the models can be dis- 
tinguished by the abundance of substructure as quanti- 
fied by the delta-deviation (100 % means 'different' and 
% means 'identical'). The values in Table || show that 
the probability distributions for the SCDM models dif- 
fer from that for the ACDM and OCDM models. On the 
other hand, the low density models with and without a 
cosmological constant are very similar, i.e., they cannot 
be discriminated. For all four models, the abundance of 
substructure depends more on the imposed limit in the 
delta-deviation than on the difference between models. 
Some degree of substructure is typical for all simulated 
clusters. 

As we have recently shown in simulations (Knebe & 
Miiller 1999), unvirialized particle groups can be identi- 
fied with recent or ongoing mergers which could be a pos- 
sible explanation for substructure. For this reason we have 
separated 'virialized' and 'unvirialized' clusters. The cor- 
responding differential distribution of the delta-deviation 
in the SCDM1 model is shown in Fig. ^. Unvirialized clus- 
ters clearly have more substructure than relaxed systems. 
Substructure is therefore a clear sign of incomplete relax- 
ation. Note that the distributions of virialized and unviri- 
alized groups are separately normalized, but the number 
of unvirialized groups is always much smaller. The same 
behaviour is also found for the other models. 



4. Comparison with Other Indicators of 
Substructure 

For testing the delta-deviation against other substructure 
indicators, we have also employed the friends-of-friends 
algorithm with a smaller linking length. Similar to the hi- 



erarchical tree (Serna & Gerbal 1996, Klypin et al. 1999) 



we determine the mass fraction of the two most massive 
subgroups to produce multiplicity statistics. To this aim, 
we build FOF clusters with half of the linking length of the 



Fig. 4. Probability distribution of A for the SCDM1 
model distinguishing virialized (solid) and unvirialized 
(dashed) particle groups (with Poisson error bars). 

original analysis. Then the total mass M of each particle 
group of the first analysis is decomposed into 



M = TOi + ?7l2 



(3) 



where mi and mi are the masses of the two most massive 
subclumps. The remaining mass m rcra consists of all the 
constituents lying in low density parts of the cluster, prob- 
ably containing mainly satellites recently accreted. We de- 
fine the multiplicity of each cluster to be 



mi + m 2 



(4) 



For M p « 1 we do not resolve substructure because the 
second most massive subgroup does not contribute to the 
total mass of the object. But for M p rj 2 both subclumps 
are of comparable mass and therefore we have a cluster 
with a possible double structure. A probable origin is a 
big merger of almost equal-mass progenitors. 

Fig. H shows the correlation of A/N and M p for the 
SCDM1 simulation. The other models look similar. The 
lines separate regions where both statistics indicate there 
to be substructure. One notes a weak correlation with a 
wide scatter and a large number of outliers. The points in 
the upper left part show high values of the delta-deviation 
for which the multiplicity test shows that the subgroups 
are not very massive. So, they are probably small satellites 
falling onto the main cluster as seen, e.g., in Fig. 0. The 
lower right part shows clusters with subclumps of nearly 
equal masses, but no velocity deviations. A possible ex- 
planation are projection effects. 

Another interesting point is the relationship of the re- 
maining mass TO rcm to the mass of the two most massive 
subclumps mi and mi. If m lcm is not negligible, we ex- 
pect a loosely bound cluster with a large portion of ac- 
creted material which may be formed recently. In Fig. |^ 
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Fig. 5. Correlation for A/N vs. multiplicity M p in the 
SCDM1 model. The lines A/N = 1.4 and M p = 1.5 show 
the adopted criteria for substructure. 

we show the correlation between m lest /mi and the mul- 
tiplicity M p for the SCDM1 simulation. The straight line 
gives the equality between mi +m2 and m rem . Above it we 
find loosely bound clusters which are obviously exceptions 
(about 13%). Below this line, the influence of recently ac- 
creted material is small. The loosely bound clusters lie 
mainly in the range where the multiplicity distribution 
indicates no significant substructure. This means that a 
steady accretion of material onto galaxy clusters tends to 
destroy substructure. Also we checked that a large accre- 
tion leads mostly to low values of A/N. 

In Fig. we show the dependence of identified sub- 
structure on the cluster mass M in the SCDM1 model. We 
find particle groups with and without substructure over 
the whole mass range, but the scatter for light clusters 
is much larger. When restricting to higher mass objects, 
the percentage of groups with substructure increases, and 
it leads to a value of about 30%, as in observed cluster 
samples (note that above we took a fixed cluster density 
which are taken from the same high mass range). 

5. Conclusions 

We quantified substructure in simulated clusters of 
galaxies using the delta-deviation statistics proposed by 
Dressier & Shectman (1988). Towards this aim, we stud- 
ied the scaling of the statistics with the multiplicity of 
the total cluster and with the multiplicity of subgroups. 
Even with an optimal choice of both parameters which 
requires a large numbers of observed redshifts, we get a 
broad distribution for the statistics. The amount of sub- 
clustering depends on the chosen criterion, and we can 
reasonably reproduce the amount of substructure found 



Fig. 6. Influence of the rest mass ™ r est on the multiplic- 
ity M p in the SCDM1 model. The straight line gives the 
equality between the masses of two biggest clumps and 
the rest. 




4x10 14 8x10 14 1.2x10 15 

M [h- ] M sol ] 

Fig. 7. Substructure dependence on mass of the galaxy 
cluster in the SCDM1 model, where the line A/N = 1.4 
gives the adopted substructure criterion. 

recently in galaxy clusters redshift surveys. The chosen 
criterion (A/N > 1.4) is similar to that used in obser- 
vations, but the scatter is high, and reasonable changes 
may lead to substructure in the range of (20 — 50)%. The 
difference between different cosmologies is significant and 
a measure of the mean matter density of the universe. 
Among the models studied, ACDM shows the smallest 
percentage, and SCDM2 the highest percentage of sub- 
clustering. We recommend the cumulative distribution of 
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the delta-deviation as a natural quantifier of substructure. 
The differences in its distribution requires large catalogues 
of galaxy clusters, each with about 50 redshift measure- 
ments. On the other hand, it is shown that subcluster- 
ing is a typical property of cluster formation in hierar- 
chical theories of structure formation. In recent hydro- 
dynamical simulations, a high percentage of substructure 
(in 4 of 10 clusters) were found in a low density ACDM 



model (Eke et al. 1998). Therefore, quantifying substruc- 



ture with a distribution function of a substructure indi- 
cator, as done here with dark matter simulations, seems 
to be a prerequisite for discriminating cosmological mod- 
els. We suspect this is a general characteristic for different 
substructure indicators. 

We compared the velocity kinematics as an indicator 
for subclusters with the subgroups found as friends-of- 
fricnds groups with half of the linking length. There is a 
weak correlation, i.e., both statistics define similar struc- 
tures, but there is large scatter. On the one hand, projec- 
tion effects influence the delta-deviation, and on the other 
hand, high values of the delta-deviation can be produced 
by small groups with separate velocity kinematics. These 
effects are often stronger than the differences between dif- 
ferent cosmological scenarios. This may be the reason for 
the large differences in the literature on the amount of 
substructure in galaxy clusters. 

It was shown that the delta-deviation is most sensitive 
to recent big mergers, and that big mergers are more typ- 
ical for the high mass clusters. Furthermore, they occur 
more often in high-density than in low-density models. 

Acknowledgements. We acknowledge the use of Couchman's 
AP3M code in this study. 
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